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We identify a new class of topologically driven phase transitions when calculating the Hall conduc¬ 
tance of two-band Chern insulators in the long-time limit after a global quench of the Hamiltonian. 

The Hall conductance is expressed as the integral of the Berry curvature in the diagonal ensemble. 

Even if the topological invariant of the wave function is conserved under unitary evolution, the Hall 
conductance as a function of the energy gap in the post-quench Hamiltonian displays a continuous 
but nonanalytic behavior, that is it has a logarithmically divergent derivative as the gap closes. The 
coefficient of this logarithmic function is the ratio of the change of Chern number in the ground state 
of the post-quench Hamiltonian to the energy gap in the initial state. This nonanalytic behavior is 
universal in two-band Chern insulators. 


Introduction.- The discovery of the quantum Hall ef¬ 
fect [1101, i.e. a quantized Hall conductance in the ground 
state which jumps from one plateau to another, inspired 
the study of topological order [1, to characterize dif¬ 
ferent topological phases outside the conventional frame¬ 
work of spontaneous symmetry breaking. Considerable 
effort has been devoted to understanding topological or¬ 
der or symmetry protected topological (SPT) order in 
the ground state. More recently, a lot of attention was 
devoted to the nature of topological order and SPT or¬ 
der for a state driven out of equilibrium, in particular for 
quantum quenches of the Hamiltonian [5l-[l3|. 


Consider an isolated system initially in the ground 
state of a Hamiltonian Hi and suddenly changing the 
Hamiltonian to Hf. The wave function follows a unitary 
time evolution, while the local observables in the long 
time limit settle to the prediction of the diagonal ensem¬ 
ble 1^ , which in some cases can be reduced to a thermal 
ensemble or a generalized Gibbs ensemble 1^, [16| . Topo¬ 
logical order or SPT order cannot be expressed as a local 
observable. Therefore, its identification in a nonequilib- 
rium state is far from trivial. In the toric code model, the 
topological entropy in the long time limit is found to be 
the same as its initial value independent from whether 
the ground states of Hi oi Hf are topologically trivial 
or not This result agrees with a universal ar¬ 

gument for gapped spin liquids [13. Similarly, for the 
Fermi gas on a honeycomb lattice which essentially sim¬ 
ulates the Haldane model, the Chern number is proved 
to be conserved under unitary evolution Blli. How¬ 
ever, in the two-dimensional topological superfluid, the 
winding number of the retarded Green’s function after a 
quench shows a str ong dependence on the post-quench 
Hamiltonian Hf [1, [J even if the winding of the An¬ 
derson pseudo spin texture is conserved [1^. Also in 
the one-dimensional case, an analysis of tunneling spec¬ 
troscopy by coupling the system to an auxiliary thermal 
bath shows that the SPT order is mostly determined 
by Hf [llj. But in topological superconductors with 


proximity-induced superconductivity, the Majorana or¬ 
der parameter [l^ or the entanglement spectrum 131 
indicate that the quenched state is topologically trivial if 
Hi and Hf are in different topological phases. 

To clarify the issue of SPT order far from equilibrium, 
we appeal to a measurable physical quantity, namely the 
Hall conductance in Ghern insulators. We first study a 
paradigmatic model, i.e. the Dirac model (^ . and then 
extend our results to a general two-band Chern insula¬ 
tor. We find that the Chern number of the unitarily 
evolving wave function is conserved and uniquely deter¬ 
mined by Hi- However, while the Hall conductance of 
the quenched state is a continuous function of the energy 
gap in Hf, the derivative of this function displays a log¬ 
arithmic divergence whenever the Chern number of the 
ground state of Hf changes. We thus identify a new class 
of topologically driven phase transitions with an exotic 
critical behavior, which is quite different from the ortho¬ 
dox one in which the Hall conductance is discontinuous 
but its derivative is zero everywhere in the phase dia¬ 
gram. The discrepancy in the SPT order obtained from 
the Chern number (based on unitary time evolution) and 
the Hall conductance is attributed to the fact that the 
latter must be calculated from the diagonal ensemble, in 
which the coherence between different eigenstates oi Hf 
in the wave function is lost in the long-time limit. In this 
experimentally relevant sense the SPT order of quenched 
states depends on Hf. 

Real-time dynamics of the Chern number.- The Hamil¬ 
tonian of a two-band Chern insulator in two dimensions 
is expressed as 


k 

where Cj: = is the fermionic operator and 

sums over a single Brillouin zone. The single-particle 
Hamiltonian can be decomposed into TL^ = df: ■ a, 
where a denotes the Pauli matrices. 
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The Dirac model is a paradigm for two-band Chern 
insulators [ 2 ^. In the Dirac model, the coefficients of 
the Pauli matrices are = {kx, ky,M — Bk^) with two 
parameters M and B, and sums over the whole mo¬ 
mentum plane. The ground state is well known to be clas¬ 
sified by the Chern number C = ^ (sgn(M) -|- sgn(i?)), 
which is quantized and changes only at the phase bound¬ 
ary M = 0 or S = 0. The Hall conductance of the ground 
state is simply the Chern number in units of e^/h. 

At the time t = 0, we suddenly change the Hamil¬ 
tonian from Hi = H{Mi,Bi) to Hf = H{Mf,Bf). 
Then the wave function evolves according to |lE'(t)) = 
= nfe®|ug(t)), where \u^{t)) is the single¬ 
particle wave function obeying 

The momentum is a good quantum number both in Hi 
and Hf. Therefore, it is natural to generalize the defini¬ 
tion of the Chern number for the time-dependent wave 
function in the following way: 


C{t) = 


i 

27r 




This real-time Chern number characterizes the topolog¬ 
ical property of the wave function |d>(t)), and can be 

where S 


reexpressed as C{t) = ^ f dS ■ x A(t 
denotes the kx-ky plane oriented in the fcz-direction and 
A{t) = {uj:{t)\V is the Berry connection. C{t) is 
determined by the poles of A{t) and must remain quan¬ 
tized at all times since locally deforming A{t) cannot 
change it. In fact, the two poles of A{t) at k = 0 and 
k = 00 have conserved residues under a unitary evo¬ 
lution [2lj, so that for arbitrary Hi and Hf we have 
C{t) = (7(0). The Chern number of the wave function 
never changes although the system is driven out of equi¬ 
librium, which agrees with the no-go theorem proved by 
D’Alessio and Rigol [l^. This result suggests that the 
SPT order of a wave function is generally conserved after 
a quench if the Hamiltonian in real space contains only 
local operators tl2|. 

Hall conductance in the diagonal ensemble.- The ob¬ 
servation that C(t) is independent ot Hf does not imply 
the absence of nonequilibrium phase transitions because 
C(t) is not a measurable physical quantity. In this pa¬ 
per, a nonequilibrium phase transition is unambiguously 
indicated by the nonanalytic behavior of observables as 
the post-quench Hamiltonian Hf changes. We choose 
the Hall conductance as the indicator of nonequilibrium 
phase transitions. Notice that in the ground state the 
Hall conductance is directly related to the Chern num¬ 
ber. 

It is well known that the Hall conductance cannot be 
expressed as the expectation value of a local operator, 
but must be written as the long-time response to an ex¬ 
ternal electric field in linear response theory. This fact 
reflects the topological nature of the Hall conductance 


and is related to the observation that in order to mea¬ 
sure the Hall conductance, one must couple the system 
to auxiliary reservoirs. However, coupling to reservoirs 
unavoidably introduces decoherence and therefore in the 
long-time limit the far-from-equilibrium system will be 
described by the diagonal ensemble and not the unitar- 
ily evolved wave function of the isolated system. This 
motivates us to pursue a definition of SPT order and 
topologically driven nonequilibrium phase transitions by 
studying the Hall conductance in the diagonal ensem¬ 
ble, which is the experimentally relevant setting. In the 
long-time limit, the off-diagonal terms of the density ma¬ 
trix in the eigenbasis ot Hf are averaged out [l^. The 
time-averaged expectation value of an operator O can be 
expressed as 

hm 1 r dt{^{t)\d\^it)) = j2\{mm\HE\d\E) 

T^oo T Jo ^ 

= Tr[Op], 

( 3 ) 


where \E) are the eigenstates oi Hf and p is diagonal 
in the basis \E) with the elements |(A|(I'(0))p. If the 
long-time limit of ('I'(t)|(7|'I'(t)) exists, it must be deter¬ 
mined by p, the so-called diagonal ensemble M- While 
this argument is based on non-degenerate eigenenergies, 
the applicability of the diagonal ensemble has also been 
shown in many integrable quantum many-body models 

min. 

We build our formalism on the diagonal ensemble with 
the density matrix written as 


E 


\a—± 


nr \u^ 

ka 


f 


ka 


( 4 ) 


where \ul ) is the eigenvector ot hI and a = ± de- 

ka' ^ k 

notes the upper and lower bands with the positive and 
negative eigenvalues respectively, is the oc¬ 

cupation number of the band a and can be expressed as 
the overlap = \{ul ImI )p, where |u7 ) is the lower- 


band eigenvector of the initial Hamiltonian "HI, which is 
in fact the initial wave function. The total occupation 
at each k is conserved to be n-I- n^_ = 1. Eq. 0 is 
obtained by averaging out the ofF-diagonal elements in 
{\Uk{t)){uj:(t)\). 

Now we calculate the Hall conductance of the diagonal 
ensemble in linear response theory 2^, i.e., we replace 
the equilibrium density matrix in linear response theory 
by the diagonal ensemble p. This replacement does not 
cause any problem in the formalism because p is time- 
independent satisfying [p^Hf] = 0. We can then express 
the Hall conductance as the current-current correlation 
in the diagonal ensemble: 


(Jh = lim / dte^ 
CJ—>-0 oLU Jq 


'Tr I 


(^P Jy: Jx{t) 


( 5 ) 
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Mj=1.01, Bi=-1.51 



Mj=1.01, Bi=1.51 



FIG. 1. (Color online) The Hall conductance Cneq as a func¬ 
tion of {Mf,Bf) at different in the Dirac model. 

[Top panel] The initial state is topologically trivial. [Bottom 
panel] The initial state is topologically nontrivial. 


where S denotes the area of the system and is conven- 

, ml ^ 

tionally set to unity. Jp = eg is the current 

operator along the /3-direction with e denoting the charge 
of the particle. Following the process for obtaining the 
celebrated TKNN number Q, we reexpress the dimen¬ 
sionless Hall conductance Cneq '■= as [21| 


a, 


neq 





dul \ 

ka \ 

dky j 



( 6 ) 

which is the integral of the weighted mixture of Berry 
curvatures in different bands of the post-quench Hamilto¬ 
nian. In the case oiHi = Hf (no quench), the occupation 
is n^_ = 1 and = 0 everywhere in the Brillouin zone, 
and Cneq is just the Chern number of the initial state as 
we expect. But tor Hi Hf, G [0,1] becomes a con¬ 

tinuous function of k so that Cneq is not quantized any 
more but can take an arbitrary value. 

It is worth comparing the real time Chern number C{t) 
in Eq. m with the dimensionless Hall conductance Cneq 
in Eq. ©• The former reflects the topology of the wave 
function, being quantized but not measurable, while the 


latter is a true observable but not quantized. They are 
both integrals of the Berry curvature, but C(t) is derived 
from the wave function while Cneq follows from the diag¬ 
onal ensemble where the coherence is lost. Decoherence 
plays a crucial role in understanding the SPT order of a 
quenched state in the long-time limit which is a nonequi¬ 
librium steady state. 

Nonanalytic behavior of Hall conductance- In the 
Dirac model, it is straightforward to determine the Hall 
conductance as 


211 


Cneq — 


dk 


(l + {Bik - Mi){Bfk - Mf)j (Bfk + Mflj 


with k + [Bi/fk — Mi/f) 


jneq is a function of {Mi, Bi^ Mf, Bf), i-e- 


H (‘‘0 

(7) 

The Hall conduc- 
the 


tance C„ 

parameters of Hi and Hf. This function satisfies the 
properties: 

Cneq {Mi,Bi, Mf , Bf ) — Cneq {B^,Mi,Bf,Mf) 

= - Cneq{-M,, -Mf,-Bf)’ 

Let us study this function as {Mf,Bf) changes, while 
{Mi,Bi) is invariant, i.e., the initial state is fixed. Due 
to Eq. ([5]), we only consider the cases Mi,Bi > 0 and 
Mi > 0,Bi < 0. As shown in Fig.[Tl Cneq is a continuous 


function of {Mf,Bf) in the whole parameter space [21 1. 
This result is surprising if we consider the fact that the 
Chern number of the ground state has a jump when¬ 
ever M or B change sign. By driving the system out of 
equilibrium, we smoothen the Hall conductance function. 
Cneq{Mf,Bf) has a similar shape at different {Mi,Bi), 
reminiscent of the function (sgn(M/) -|-sgn(i3/))/2, i.e., 
the Chern number in the ground-state wave function 
of the post-quench Hamiltonian Hf. As Mf,Bf 0 
{Mf,Bf <C 0), Cneq takes a positive (negative) value, 
while Cneq is close to zero as Mf and Bf have different 
signs. Even if the initial state is topologically trivial (see 
Fig.lll the top panel), the Hall conductance is finite as Hf 
is in the nontrivial regime, but it cannot reach the quan¬ 
tized values ±e^//i. When the initial state is nontrivial 
(see Fig. [H the bottom panel), the Hall conductance is 
suppressed as Hf deviates from Hi, and can even change 
the sign as Mf and Bf both change their signs. 

While Cneq{Mf,Bf) is continuous, the key finding is 
that whenever the post-quench Hamiltonian crosses the 
boundary at Mf = 0 {Bf = 0), the derivative of the Hall 
conductance ) diverges to -boo in a logarith¬ 


mic way 2l|: 


dC, 


neq 


1™ 

Mf-iO dMf 

7)Cneq 
Inn ——- 
Bf^O oBf 


-1 

-1 


ln|M/|, 


( 9 ) 


2LB,; 


■ In [Hf 
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As Mf —>• 0, as a function of (ln|M/|) asymp¬ 

totically approaches a straight line with the slope 

— l/(2|Mi|), which is independent of Bf,Bi and the side 

dC 

from which Mf goes to zero. As Bf —?► 0, has a 

similar divergence since Cneq is invariant under the ex¬ 
change oi Mf and Bf according to Eq. ©■ We identify a 
nonequilibrium phase transition when the Chern number 
in the ground state oi Hf changes. The critical behav¬ 
ior of this phase transition is exotic, compared to that 
of ground-state phase transitions in which the Hall con¬ 
ductance has a zero derivative everywhere but displays a 
discontinuity at the phase boundary. 

This phase transition reveals different nonequilibrium 
phases which share the common symmetries of the Dirac 
model. Apparently, the broken symmetry picture does 
not account for this transition, which must be topolog¬ 
ically driven. Interestingly, the topological invariant of 
the wave function C{t) is independent of {Mf,Bf), and 
then fails to characterize different phases in this nonequi¬ 
librium phase transition. One can assign the Chern num¬ 
ber C{Hf) of the ground-state wave function oi Hf to 
each nonequilibrium phase to distinguish them. We will 
see that the change of C{Hf) determines the character 
of this nonequilibrium phase transition in a general two- 
band Chern insulator. 

Now let us consider a general two-band Chern insu¬ 
lator in two dimensions with the Hamiltonian given by 
Eq. ([T]). The coefficient vector 

different from model to model. But the nonanalytic be¬ 
havior of Hall conductance is insensitive to the change of 
d^. Instead, it depends only upon the lowest-order ex¬ 
pansion of dj: at the momentum q where the energy gap 
closes {dq = 0) at a phase boundary. In a generic model, 
two components of d^ must be zero. Let us suppose them 
to be dig- and d2q without loss of generality. The energy 
gap is dq- = Ida^l when the system is close to the phase 
boundary, d^^ is a free parameter in the Hamiltonian 
(the gap parameter), which is denoted by m. Note that 
TO = M in the Dirac model. 

Suppose that the system is initially in a ground state 
with the gap parameter to = TOi, before we suddenly 
change to in the Hamiltonian from rrii to mf. We mea¬ 
sure the Hall conductance in the long time limit. The 
Hall conductance Cneq is a function oi mf, while we fix 
TOi to be nonzero. The function Cneq(jnf) is continu¬ 
ous but nonanalytic at mf =0, where the gap of the 
post-quench Hamiltonian Hf closes. The derivative of 
Cneqijnf) satisfies [21 1 


lim 

dmf 


lim Cijrif) — lim C{mf) 

my—^0+ 


2|TOi| 


ln|TO/|, 


( 10 ) 

where C{mf) denotes the Chern number in the ground- 
state wave function oi Hf. The derivative of the Hall 
conductance with respect to the energy gap in dfy is log¬ 


arithmically divergent as the gap closes. And the co¬ 
efficient of this logarithmic function is the ratio of the 
change of Chern number in the ground state oi Hf to 
the energy gap in the initial state. Eq. m relates the 
nonequilibrium phase transition in quenched states to the 
topological phase transition in ground states, indicating 
that this nonequilibrium phase transition is in fact topo¬ 
logically driven. Eq. for the Dirac model is a special 
case of Eq. cni) as the change of Chern number is —1. 

Conclusions.- In summary, we find a new class of topo¬ 
logically driven phase transitions in quenched states of 
two-band Chern insulators, which are characterized by 
the Hall conductance as a continuous function of the 
energy gap in the post-quench Hamiltonian Hf with a 
logarithmically divergent derivative. The asymptotic be¬ 
havior of the Hall conductance function is determined by 
the ratio of the change of Chern number in the ground 
state oi Hf to the energy gap in the initial state, which 
is universal in two-band Chern insulators. We obtain 
the Hall conductance by applying linear response the¬ 
ory in the diagonal ensemble of the system, which is the 
physically correct description of the long-time limit in 
a far-from-equilibrium quench setup. The topological in¬ 
variant of the real-time wave function fails to predict this 
phase transition, which can only be correctly identified 
in the diagonal ensemble where decoherence effects are 
taken into account. Our finding indicates the possibil¬ 
ity of exotic topological phase transitions in systems far 
from equilibrium. 


Finally, we discuss the conditions for observing this 
phase transition in experiments. The nonequilibrium dis¬ 
tribution of particles is responsible for the logarithmi¬ 
cally divergent derivative of the Hall conductance. Ul¬ 
tracold atomic gases are known to be well isolated from 
the environment and suitable for studying the quench 
dynamics of many-body quantum systems [25l| . The Hal¬ 
dane model [2^ was recently realized with cold atoms 
in an optical lattice The Haldane model is a 

two-band Chern insulator, in which the quenched-state 
Hall conductance displays the nonanalytic behavior in 
Eq. (flUl) [2^ . The measurement of conductances in cold 
atoms is difficult^but a two-terminal setup was imple¬ 
mented recently 2^, 3^. We expect that our prediction 
can be checked in a four-terminal setup made of cold 
atoms simulating the Haldane model. 
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SUPPLEMENTARY MATERIAL 


The real-time Chern number C{t) in the Dirac model 


We express the real-time Chern number as 

C{t) = ^ fdS-(vf:xA{t)), ( 11 ) 


where A{t) = {u^{t)\Vj^\uj^{t)) is the Berry connection 
with \uj:{t)) = {(pi, p 2 )'^ denoting the single-particle wave 
function. In the Dirac model, it is straightforward to 
calculate the wave function and obtain 




1 


-:-^-= [ - COs{dU){d\ - d 


-I- i sin(d^f)- 


3fc 


(d^ -P 

\ k 3kJ 


( 12 ) 


and 


(p2{t) =- 


2diJ2dUdi-d^A 

ky k^ k 3k^ 


[2di cos{dlt) 


2i sin(dit)(dp — dt 


where = kx + iky, d 


k ^ ^ k 

?/f 


3k 


■ 4 )]> 


(13) 


= (d‘C,d‘C,d‘C) is the co- 

^ U- ’ 2k ’ 3fc ' 

d 


efficient vector in the initial and post-quench Hamilto- 

Hi ~il f 

nians, respectively, and d- is the length of d- . We 
divide the Berry connection into A{t) = Ai{t) -|- A 2 {t) 
with Ac = (/)* Vj:d>a- Noticing that pi is a function of 


k = ijk‘^ + ky, we immediately know that x 

must be zero, so that Ai does not contribute to C{t). 
We again divide A 2 into the irrelevant term A 2 a with a 
zero curl and the relevant term A 2 b with its imaginary 
part written as 


Im[A2b] 


(dt 

cos'^ {dlt) +sin^ (dlt)—— 

k k 


— kyX - 


kxV 


2dt{di - dU) ’ 

k 3k' 



(14) 


where x and y denote the unit vectors in the momentum 
plane. 

Now we reexpress the Chern number by the vector field 
Im[A 2 b] as 

C{t) = ^JdS- (vg X Im[A26(<)]) . (15) 

Im[A 2 f,] is a vortex field with two poles at zero and inhn- 
ity, respectively. Applying the Kelvin-Stokes theorem in 


an annulus with inner radius r and outer radius R, and 
then taking the limit r —> 0 and i? —>■ 00 , we obtain 

—27rC'(t) = lim f dS ■ (vr X Im[A 2 b]') 

= ( lim ® — lim ® I Im[A 2 bl • dk 

\^^°°Ik=R ^^'^Jk=rj ( 16 ) 

= lim (27ri?|Im[A26]|fc=/?) 

- lim (27rr|Im[A2b]|fe=r ) , 

where |Im[A 2 b]|fe=H denotes the length of the vector 
Im[A 2 f,] at the circle of radius {k = R). The first limit 
evaluates (7r(l —sgn(i3i))), while the second limit evalu¬ 
ates ( 71(1 -I- sgn(Mi))), being both time-independent. In 
other words, the residues of Im[A 2 b] at zero and infin¬ 
ity are both time-invariant, which leads to a conserved 
Chern number: 

Cit) = ^ (sgn(M,) -f sgn(B,)). (17) 


The Hall conductance of quenched states 


In this section, we first show how to express the Hall 
conductance of quenched states as the integral of the 
Berry curvature. Our derivation is a straightforward ex¬ 
tension of the work by Thouless et al. [l|. We then ex¬ 
press the Hall conductance by using the coefficient vec¬ 
tors in two-band Chern insulators. 

In linear response theory, the Hall conductance is writ¬ 
ten as 


an = lim — 

cj—^■O OLO 


^ dtd^^-^d\T^[p[jy,Jx{t)]), (18) 


where 77 is an infinitesimal number corresponding to the 
adiabatic switch-on of an external electric field, and w is 
the frequency of the electric field with the limit w —>■ 0 
corresponding to the dc-conductance. The diagonal en¬ 
semble is known to be p = Jlfc® (Za ’ 

which is a product state. Due to the conversation law 
Sa '^ka ~ state of the system is limited in a sub¬ 

space of the Fock space in which the empty or doubly- 
occupied states at each momentum are excluded. We can 
then reexpress an in the first-quantization language as 


1 r°° 

an = lim — / 
cj—^■O oLO Jq 


X } nr { u\. 

kOL \ fco; 


Jim) 


kOL 


(19) 


where the momentum-resolved current operator is J- = 

ml 

e ——. Since we are interested in the dc Hall conductance 

okp 
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which is a real number, we take the real part of <jh and 
obtain 


RecT// = 


(JH + (T 


H 


—le 


E 

k,a,^ 


^kot 


(ej — e 

k kn/ 


koL 


koL 

dnl 
_^ 

dky 


kp) 


fe/3 


fe/3 


dnl 
_^ 

dk^ 


koL 


-H.c. 


( 20 ) 


where denotes the eigenvalue of T-LI. We make use of 
the relation T-it = et \ui ){ul \ and finally obtain 

k ^ kcn kcn 


a 


neq 


ReCT// 

jh 


% 

27r 


E 




dui \ 

koi \ 

dky j 



In a two-band Chern insulator with the Hamiltonian 
= dg • (T, the Berry curvatures in different bands are 
opposite to each other. By using this and the conserva¬ 
tion law -I- ng_ = 1, we reexpress Eq. m as 


case at My = 0 when By is fixed to be nonzero, since 
Cneq{Mf, Bf) IS invariant under the exchange of My and 
Bf- 

In the Dirac model, both C and cos 9 are rotationally 
invariant in the kx-ky plane. We can then carry out the 
azimuthal integration in Eq. (I22|) . By making a substi¬ 
tution k = k^, we express the Hall conductance as 


pOO 

Cneq = / dk COS 9 ■ C, 

Jo 

where the Berry curvature is 


1 B fk + Mf 


(25) 


(26) 


and the occupation factor is 


didl 

k k 


with = ^Jk + {Bi/fk — Mi/fY. At My ^ 0, we can 
express the derivative of Cne.q as 


r^~d{cos9-C) 
dMf Jo dMf 


(28) 


C'neq — 


dk^ COS 9 ■ C, 


( 22 ) 


where C denotes the Berry curvature in the lower-band 
of the post-quench Hamiltonian H y and can be expressed 
as 


ddt ddl 

_fc ^ _fc 

dkx dky 

Mdir 

and cos 9 is the occupation factor defined as 




(23) 


cosd :=(2ng_ — 1) 

=(<;]. <?e)/(<4) 


(24) 


with 9 denoting the angle between dt and dl. dt and 

k k k 

dt are the coefficients of the Pauli matrices in the initial 

k 

and post-quench Hamiltonians, respectively, and dt and 
dl are the length of dt and dl, respectively. 


Continuity and nonanalyticity of the Hall 
condnctance in the Dirac model 

In this section, we show how to prove the continuity 
of Cneq (My, By) and the logarithmic divergence of its 
derivative at the phase boundary. We only prove the 


A straightforward observation is that both cosd(fc) and 
C{k) are smooth functions for k G (0, oo). How¬ 
ever, they do not uniformly converge to (cosd)M/=o 
or Cmj=o as Mf —>■ 0. The unique singularity is 
fc = 0, at which we have lim^_,^Q limM^-s-o cos d = 0 
but liniM/^o linr^_^Q cosd = sgn(Mi)sgn(My). And 

C (L = (Jj = sgn(My)/(4M|) is divergent as Mf —>• 0. 

We divide the integral into two parts: dk = dk + 

dk with 77 > 0 a number that can be arbitrarily small, 
xle second integral is a smooth function oi Mf, which 
can be proved by studying the asymptotic behavior of 
(cosd -C) in the limit fc —> 00 , or more precisely, by mak¬ 
ing a substitution /c —>■ 1/fc in the integral. In fact, k = 00 
is a true singularity at the boundary By = 0, where k = 0 
is a regular point, since cos 9 and C are invariant under 
the substitution k l/k and Mjyy B^jf. If there is 
any nonanalytic behavior in the function Cneq{Mf), it 
must come from the first integral denoted by next. 
Interestingly, we can choose an arbitrarily small 77 so that 
dk in cosd converges to a constant dk = \Mi\. We then 
obtain 

fv (k + {Bik-Mi){Bfk-Mf)){Bfk + Mf) 

Cneo = / dkk^ - - - 9 - 

•^0 4|M,| (A: + (Byfc-My)2j 

(29) 

The calculation of this integral is straightforward since 
the integrand is rational. 
























is lengthy, but it is an elementary function. F{r]) is a 
smooth function of M/, while F{0) is expressed as 


We express the result as C^eq = F{r]) — F{0) with 
F denoting the original function. The expression of F 

I 


F(0) = 


+ Bf)Mf - ^ - ^BfMf + 1^ y/1- 4:BfMj - SBjMj + 6BfMf - 1 


S\M,\Bj [ ^l-ABfMf 


Mj 


(30) 


+ B, lnM| 


2B, - 2B 


We are interested in F{0) as a function of M/ in the 
neighborhood of the phase boundary Mf = 0. We notice 
that — ABfMf can be expanded at Mf = 0 into 


V^l - 4%M/ =1 - 2BfMf - 2BjMj - 4:BjMj 
- lOBjMf + 0{M^). 


(31) 


We substitute this expression into Eq. o and obtain 


F(0) 


2B,-2Bf- B,\n{2Bj) Mf . 

8\M,\B} 4|M,| f 

+ 0{Mf) In (1 + 0{Mf)) + 0{Mf) In \Mf\. 

(32) 



FIG. 2. dCneq/dMf as a function of (ln|M/|) at different 
(Mi, Bi, Bf) in the Dirac model. Note the curves at Mi = 1, 
in which we simultaneously plot the data at Mf ^ 0'^ and at 
Mf —>■ 0“, which are in fact undistinguishable at small \Mf\. 


The first term is independent ot Mf. The second term is a 
continuous function oi Mf, but its derivative with respect 
to Mf is divergent as Mf 0. All the other terms are 
continuous functions ot Mf, and their derivatives with 
respect to Mf are finite at M/ = 0. The asymptotic 
behavior of dCneq/dMf is uniquely determined by the 
second term. The function CneqiMf) then asymptoti¬ 
cally approaches (—M/In |M/|/(2|Mi|)-b const.) in the 
limit Mf ^ 0. This immediately leads to our results 
that CneqiMf) is continuous 0 and dCneq/dMf is loga¬ 
rithmically divergent as 


lim 

Mf^ 


dCn 


0 dMf 


-1 

2\M~\ 


ln|M/|. 


(33) 


Furthermore, we calculate the Hall conductance by nu¬ 
merically integrate Eq. (1^ . We plot as a func¬ 

tion of ln|M/| in Fig. [21 In the limit Mf —>■ 0, the 
curves asymptotically approach straight lines with the 
slope —1/21 Mi I, which is independent of i?i, Bf and the 
side from which Mf goes to zero. The numerical result 
coincides well with our analysis. 


Universal nonanalytic behavior of the Hall 
conductance function in two-band Chern insulators 

Let us consider a general two-band Chern insulator in 
two dimensions with the Hamiltonian expressed as 

^ = ( 34 ) 

k 

where the single-particle Hamiltonian can be decomposed 
into ■ S with cj denoting the Pauli matrices. 

Examples include the Dirac model, the Haldane model 31 
or the Kitaev honeycomb model in the fermionic basis 4, 
Q. The coefficient vector dg = {d-ij:, d^j^, d^j:) is different 
from model to model. But the nonanalytic behavior of 
Hall conductance is insensitive to the change of dg, but 

depends only upon the lowest-order expansion of d^ at 
the singularities of the Berry curvature. 

Let us first show how the Chern number of the ground- 
state wave function is related to the expansion of djj. The 
Chern number is expressed by the Berry connection as 

/d5-(vgxImA) (35) 






















9 


with 


Eq. (1551) in the limit ry —>■ 0, which evaluates 


Tm/l’ = k^lk 

- d,^) • 


(36) 


C-l?! 


- (1 + sgn(m)) sgn [aixa 2 y - a 2 xaiy) ■ (41) 


The Chern number must be zero when A has no singu¬ 
larity in the Brillouin zone. According to Kelvin-Stokes 
theorem, the Chern number can be expressed as the line 
integral of A over the boundaries of the infinitesimal 
neighborhoods of singularities. Suppose that A has a set 
of singularities qi,q 2 ,--- ,qN in a single Brillouin zone. 
The Chern number can be expressed as 

N 

= (37) 

i=i 

with 

(X ImA • dk, (38) 

2tT 

where dBj^{qj) denotes the boundary of a circle of radius 
centered at qj , and the integral is along the anticlock¬ 
wise direction. Here we do not consider the singularity 
at infinity, since the Brillouin zone is finite in a generic 
model. 

In general, a singularity of A is a point q at which 
(dq-— ^ 3 ^ = 0) and then di^ = d 2 q = 0. In a generic 
model, dg- = jdag-j is the energy gap when the system is 
close to the phase boundary, d^^ is a free parameter in 
the Hamiltonian, which is denoted by m next. Note that 
TO = M in the Dirac model, m is zero if and only if the 
energy gap closes accompanied by a change of the Chern 
number. 

The Berry connection can be reexpressed as 


ImA = 


dk ^3fc ^ 

2d,- ) 


^ifc^fc^2fc d2k^kd^^ \ 


(39) 


Since dg* and d^^ are finite at to 0, we can replace 
value at fc = (f, which is (I -|- sgn{m))/2 

with sgn(TO) denoting the sign of to. This replacement 
will not change the integral in Eq. (l38ll in the limit 77 —?► 0. 
The value of d^j: a,t k ^ q has nothing to do with the 
Chern number. 

From Eq. (1551) . we know that the Chern number de¬ 
pends only upon d-j^ around the singularities of A. We 
then expand d^^ and d^j: into power series of Afc = k — q. 
Without loss of generality, we have 

d;^g =aixAkj: + aiyAky + 0{Ak‘^), 

d^j: =a2xAkx + a2yAky + OiAk"^), (40) 

dgg =TO -I- 0{Ak). 

It is straight forward to prove that the higher-order terms 
in this expansion do not contribute to the integral in 


It is worth mentioning that the three components of d^ 
are on an equal footing. Depending on the basis that 
is chosen, the components of d^: could be exchanged in 
some models. 

Notice that, in Eq. (HOI) . the coefficients oia,, aiy, a 2 x 
and a 2 y are g-dependent. While to at different qj may 
represent different parameters in the Hamiltonian, i.e. 
the gap parameters at different phase boundaries. An 
example is the Haldane model [3|. In a single Brillouin 
zone, A has two singularities. And the energy gap closes 
at one of them as the system is at some phase boundary, 
but closes at the other singularity as the system is at 
the different phase boundary. On the other hand, if the 
system has some symmetries so that at a specific phase 
boundary the gap closes simultaneously at several g}, to 
at these g) must be the same parameter. 

Now let us discuss the Hall conductance of quenched 
states when the parameters in the initial and post-quench 
Hamiltonians are both nearby a specific phase boundary 
where the gap parameter is denoted by to. Suppose that 
the system is initially in a ground state with the gap pa¬ 
rameter TO = TOi, before we suddenly change to in the 
Hamiltonian from rrii to to/. We then measure the Hall 
conductance in the long time limit. The Hall conduc¬ 
tance Cneq N a function of to/, while we fix rrii to be 
nonzero. 

Noting = d^{mi) and dZ = dj^{mf ), we express the 
Hall conductance as 


Cn 


1 

47r 


( ddd^f) ddd^f) \ 

- y dkx dky j 

{dk{mf)Y 
^ dki.mi) ■ dj:{mfA 
dkirrii) \ ’ 


4(™/) 


(42) 


where the integral is over a single Brillouin zone. In 
a generic model, the components of d^ are all analytic 
functions of k. According to Eq. (I55|) . Cneq{mf) is non- 
analytic only if d^{mf) in the denominator of the integral 
vanishes at some fc, i.e., the singularities of the Berry cur¬ 
vature. This is the case at to/ = 0 when the gap of the 
post-quench Hamiltonian closes at some singularities of 
the Berry connection A. Without loss of generality, we 
suppose that these singularities are gi, g 2 , • • • , gV' with 
N' < N. The nonanalyticity of Cneqinif) comes from the 
integral over the neighborhoods of gi,g 2 ,-- - ,qN>- We 
then divide Cneq into the analytic part and the nonan- 
alytic part as we did in the Dirac model. The latter is 
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written as 

N' 

= (43) 

i=i 

with 


niqj) 

^neq 


1 

47r 



dl? 


( ddi^jmf) dd^{mf) \ 

y dk^ dky j 


d^{mi) • dj:{mf) 
d^{mi) 


dki^f) 


(44) 


where By {qj ) is a circle of radius centered at qj with 
r] a positive number that can be arbitrarily small. 

In the neighborhood of the singularity q, we can ex¬ 
pand the components of d^ into power series. Let us first 
consider the lowest-order term given by Eq. (HUl) . We 
substitute Eq. (HHl) into Eq (Hill . We replace d^{mi) by 
its value at k = q, that is d^{mi) = \mi\. This replace¬ 
ment will not change the nonanalytic behavior of C^q 
since rrii is nonzero. While the denominator of the inte¬ 
grand becomes 

^ / 2 ^ 

— j ^ ^ ^ i^ajx^kx ~\~ ^jy^^y) 

V j=i 

(45) 

We change the coordinate system so that the function 
(^di^{mf)^ has rotational symmetry around q. In the 
new coordinate system we have 

2 

'^{ajx^kx + ajyAky)'^ = Ak'"^. (46) 

i=i 

This transformation is always possible. Otherwise, the 
coefficients before A/c^ and Aky have different signs, 
which contradicts the proposition that q is an isolated 
singularity. In the new coordinate system, we carry out 
the azimuthal integration and obtain 

/^(q) ^2x^ly) 

rriimf + Ak'^ (47) 

^m'j + AA:'2^ 

In the numerator of the integrand, only the 2nd-order 
term Afc'^ has a contribution to the nonanalyticity of 
C^q{mf). It is trivial to find the original function of 
this integral, whose value is an analytic function of mj 
at Ak'^ = Tj but a nonanalytic one at Ak'^ = 0. This 
coincides with our expectation that the nonanalytic be¬ 
havior of Cii%{mf) should be independent of the choice 


X r d{Ak''^) 

Jo 



of Tj. The nonanalytic part of C^qijnf) is 

c®, ~ (48) 


Eirst, C^q is a continuous function of nif, and then the 
Hall conductance Cneyi^nf) must be continuous. Second, 
the derivative of C^q with respect to to/ is logarithmi¬ 
cally divergent in the limit to/ —>■ 0, i.e.. 


,. dCneq Sgn(ni3;a2y ^2a:^ly) , i i /' 4 n^ 

lim —-^ , - ^InTOf. (49) 

m /—>0 dnif 2|TOq 

Comparing Eq. (HT|) with Eq. (H^ . we immediately find 
that the (f-dependent coefficient in is equal to the 

change of C^^{mf) at the phase boundary m/ = 0. 

is the sum of c!i^q at the singularities qi,q 2 ,--- ,qN>j 
while the Chern number C is the sum of at all the 
singularities of A. But at j > N' does not change 
at TO/ = 0, since the corresponding gap parameter is 
different from to/. We hnally obtain 


lim 

m/—)-0 


dC, 


neq 


druf 


lim C{mf) — lim C{mf) 

mf—^0~ rrif—^Q'^ 


2\mi\ 


ln|TO/|, 

(50) 


which is the central result of this paper. 

Eq. (1501) is obtained by considering only the lowest- 
order term in the expansion of dj^. Next we prove that 
the higher-order terms do not change the continuity of 
Cneq Or the asymptotic behavior of dCneq/drrif in the 
limit TO/ —> 0. This is true if the higher-order terms 
do not change the continuity of cj^q or the asymptotic 
behavior of dci^q / dm / at an arbitrary singularity. 

A linear term like {a^xAkx + a^iyAky) is not allowed 
in the expansion of d^^: in Eq. (HO)) . Otherwise, d^ is 
not the energy gap, or the minimum point of dj: is not 
at k = q, but changes with to, which contradicts our 
proposition. In a generic model like the Dirac model, 
the Haldane model or the Kitaev honeycomb model, the 
minimum point of dj^ is determined by the symmetry of 
the model and then keeps invariant as the system is in 
the vicinity of the phase boundary. 

Let us add the 2nd-order term into i.e., 

ipOxAkx + b^yAky + b'irnAkxAky'j without loss of gener¬ 
ality. The denominator in the integrand of C^q becomes 


(dki^f))'^ = 


^ ( {(^jxAkx 


XjyAky) 


i=i 


-I- 2to/ [bsxAkl + bsyAky + bsmAkxAky) 


+ 0 {Ak^) 


( 51 ) 
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Cnfq is an integral over the infinitesimal neighborhood of 
q, where the 4th-order term O(Afc^) is much smaller than 
the 2nd-order term and can be neglected. At the same 
time, the additional 2nd-order term that is proportional 
to nif has no contribution to the asymptotic behavior of 
C^q and dC^q/druf in the limit m/ —?► 0. Therefore, 
the effective denominator is the same as Eq. (|4^. The 
numerator of the integrand becomes 


y dk:c dky j 




{dj^{mi) ■ d^{mf)j 


— {Ojlx(^2y (^2xO>ly') 


rriim 


2 

/ 


2 

+ m/ E {^(Xjx^^kx ^jy^^^y'} 

3 = 1 


+ rrij (b^x^k^ + b^y/S.k^ + b^m^kx^ky') + 0{/\k'^) . 

(52) 


The 4th-order term O(Afc^) can be neglected in the limit 
?7 —>■ 0. This can be easily verified by adding Afc'^ in the 
numerator of the integrand in Eq. (1471) and checking the 
output. The additional 2nd-order term that is propor¬ 
tional to m'j leads to a correction of C^q ~ mj In |m/|, 
which does not change the asymptotic behavior of C^q 
and dC^q/dmf in the limit m/ —>■ 0. In the power se¬ 
ries of dgg, any term in order higher than 2 leads to a 
correction to numerator or denominator of the integrand 
which is at least in the 3rd order of Afc and can then 
be neglected in the limit 77 —>■ 0. Therefore, the higher- 


order terms in do not affect the asymptotic behavior 
of dC^q / dmf or the continuity of C^q. 

Similarly, we can prove that the higher-order terms in 
d^^ or c? 2 fc have no contribution. In fact, the terms in 
order higher than 1 lead to a correction of 0{Ak^) in 
the denominator. The terms in order higher than 3 also 
lead to a correction of O(Afc^) in the numerator, which 
can be neglected. The 2nd- and 3rd-order terms in d^^ 
or d^^ generate a linear term and a 2nd-order term that 
is proportional to m'j in the numerator. The latter does 

not contribute to the asymptotic behavior of dC^q / dm / 
due to the similar reason mentioned above. While the 
linear term in the numerator is an odd function of Akx 
or Aky, and then has no contribution to the integral since 
both the denominator and the integration boundary have 
rotational symmetry with respect to the singularity. 
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